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Puzzling discrepancies between recent pion interferometry data on Au + Au reactions at s/s = 130 
and 200 AGeV from RHIC and predictions based on ideal hydrodynamics are analyzed in terms of 
covariant parton transport theory. The discrepancies of out and longitudinal radii are significantly 
reduced when the finite opacity of the gluon plasma is taken into account. 
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Introduction. Decoupling, or freezeout of final state 
interactions, is a process with unique sensitivity to the 
space time evolution of hadronic interactions. It reflects 
the interplay between the decreasing opacity of the sys¬ 
tem at late times and correlations induced by collective 
expansion driven by high scattering rates at early times. 

In heavy-ion physics, information about the spacetime 
decoupling geometry can be obtained via identical boson 
(Hanbury-Brown and Twiss (HBT)) interferometry. Re¬ 
cent two pion correlation data in Au+Au at ^fs = 130 
JT], U and 200 AGeV/nucleon at RHIC seem at first 
sight to indicate a sudden freezeout that is difficult to 
reconcile with the strong collective dynamics implied by 
the substantial elliptic flow V 2 ~ 0.1 also observed [|] |||. 

Ideal Euler hydrodynamics provides one of the power¬ 
ful covariant approaches to predict the collective flow pat¬ 
tern as well as possible freezeout hypersurfaces in heavy- 
ion collisions. However, this approach predicts an ’’out” 
radius, R ou t as defined below, significantly larger than 
the ’’side” radius, -R S jd e [[| 0- On the other hand, it is 
possible that the failure of non-dissipative hydrodynam¬ 
ics to correctly describe the delicate space-time decou¬ 
pling geometry is due to the neglect of deviations from 
local equilibrium throughout the evolution. HBT pre¬ 
dictions with hydrodynamics are based on an additional 
ad hoc (Cooper-Frye) postulate that freezeout occurs on 
some ’’thin” three-dimensional hypersurface (typically, 
an isotherm). Only a detailed covariant transport the¬ 
ory can assess the theoretical error introduced by this 
postulate. 

The simplest Lorentz covariant dynamical framework 
that can predict freezeout self-consistently is transport 
theory^, [)| |l(], [dJ. In this approach the interaction rate 
is controlled by microscopic differential cross sections, da. 
As the system expands and rarefies, the scattering rate 
decreases until the particles stop interacting. In || it was 
emphasized that finite cross sections are needed to ac¬ 
count for the saturation of differential elliptic flow v -2 (p±) 
observed for p± > 2 GeV in Au+Au at RHIC0, |. 

The influence of final state dissipation on HBT was 
studied in a hybrid hydrodynamic/transport model in 
Ref. | p2| . In that work, hydrodynamical evolution was 
followed only up to a hadronization isotherm T(a+) = T c . 


Subsequently, the decoupling of the hadron gas was com¬ 
puted via the UrQMD hadronic transport model. How¬ 
ever, the predicted R ou t /Rside > 1 increasing with trans¬ 
verse momentum still fails to account for the observed 
decreasing R ou t/Rside < 1- Similar hadronic transport 
results were reported in jT|] . This suggests that possible 
deviations from the local equilibrium assumption prior 
to hadronization should also be explored. 

Recently, a combined parton/hadron transport theory 
approach was proposed in Q. The results suggest that 
the HBT radii are indeed sensitive to the parton cross 
section during the dense partonic phase of the reaction. 
Unfortunately, all the transport calculations above left 
open the delicate question of Lorentz covariance of the 
numerical solutions. In any case, we note that no trans¬ 
port or hydrodynamic calculation has as yet been able 
to reproduce the phenomenological decoupling source 
parametrizations fitted to that data in [fj~5| . 

In this letter, we concentrate exclusively on the par- 
tonic transport phase to isolate more clearly the influ¬ 
ence of dissipative partonic processes on the decoupling 
geometry and study in detail the question of covariance. 
Preliminary results were reported in p+| |. We utilize the 
MPC numerical technique^^, and compute the covari¬ 
ant freezeout distributions for a wide range of RHIC ini¬ 
tial conditions as a function of the partonic opacity. 

Two-particle HBT interferometry. For a chaotic bo¬ 
son freeze-out source p(x , t), the two-particle momentum 
correlation function C{p\,p 2 ) is given by the space-time 
Fourier transform of the source |u| [li], ^o). Although 
the relation cannot be inverted (due to loss of phase in¬ 
formation and the on mass-shell constraint), HBT mea¬ 
surements provide a unique test of freezeout distributions 
and dynamical scenarios in heavy-ion collisions. 

Conventionally, the two-particle correlation function is 
expressed in terms of the relative momentum q^ = p^ —ptf 
and average pair momentum K M = (p^ + P 2 V 2 . The 
’out-side-long’ variables, qo , qs , and qL , are then defined 
via qhGMS = (lO: QSt Ql) in the Longitudinal Center of 
Mass System (LCMS) reference frame where Kg CMS = 
(K°,K_ l, 0,0). We denote here the LCMS spacetime co¬ 
ordinates by x^ LCMS = (t,xo,xs,x L )- 

Experimentally, the measured correlation function (af- 
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ter correcting for Coulomb distortions) is fitted with a 
Gaussian, which for central collisions and at midrapidity 
is constrained by symmetry to the form 


C(q, K) = 1 + X(K) exp 


E ii R iW 

i=0,S,L 


(1) 


Here Ro , Rs and Rl are the ’out’, ’side’, and ’long’ HBT 
radii. For a perfectly Gaussian source, the correlation 
function is Gaussian and J2l| 

Ro(K) = {Axq)k +v±(AP) k - 2v±(Ax 0 At) K 
R 2 S (K) = (A x 2 s ) k , and R 2 L (K) = (Ax 2 L ) K , (2) 

where v± = K±/K°. Thus Rs and Rl have simple ge¬ 
ometric interpretation as the ’side’ and ’long’ widths of 
the distribution function, while Ro is a mixture of the 
’out’ width, time spread, and the xo — t correlation. 

Covariant parton transport theory. We consider here, 
as in Refs. §10- the simplest but nonlinear form of 
Lorentz-covariant Boltzmann transport theory in which 
the on-shell phase space density f(x,p ), evolves with an 
elastic 2 —> 2 rate as 

Pidfj.fi = S(x,pi) + y JJJiyhh ~ / 1 / 2 ) \Mi2^34\ 2 
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X S 4 (pi+p 2 -p 3 -P4 ) • (3) 

Here |A4| 2 is the polarization averaged scattering ma¬ 
trix element squared, the integrals are shorthands for 
Si = f d 3 Pi/[(2n) 3 Ei], while fj = (2tt) 3 f(x,pj). The 
initial conditions are specified by the source function 
S(x,p). For our applications below, we neglect quark 
degrees of freedom and interpret /(x,p) as describing an 
ultrarelativistic massless gluon gas (8 colors, 2 helicities). 

Eq. (j|) can be extended to include inelastic ma¬ 
trix elements, such as gg <-> ggg , and proper Bose or 
Fermi statistics can be introduced as well. However, at 
present there is no practical algorithm to compute accu¬ 
rate numerical solutions to such transport equations on 
the workstations available to us. Therefore, the present 
study is limited to the classical case with elastic 2 —> 2 
interactions. 

The elastic gluon scattering matrix elements in dense 
parton systems were modeled with the isotropic form 
da e i/dt = ao(s)Js, as justified by our previous study ||. 
We showed in 0] that the covariant transport solutions 
do not depend explicitly on the differential cross section 
but only on the transport opacity 

X =——( n ) ~ a tr{ dzp (x. 0 + zh.T =-)) , ( 4 ) 

Tel J V c) 

where Otr (s) = f da e i sin 2 d cm is the transport cross sec¬ 
tion (in our case, at r = 2<ro/3), and (n) is the average 
number of scatterings per parton. For a fixed nuclear 
geometry, a given transport opacity \ represents a whole 


class of initial conditions and partonic matrix elements, 
as demonstrated by the approximate proportionality || 
X oc a tr dN g (T 0 )/dg. 

We solved Eq. numerically via the MPC parton 
cascade algorithm [pd]. MPC utilizes the particle subdi¬ 
vision technique]^], which is essential to eliminate nu¬ 
merical artifacts caused by frame-dependent collision or¬ 
dering and acausal (superluminal) propagation due to 
action at a distanceJT(J. For initial partonic densities 
expected at RHIC, the severe violation of Lorentz co- 
variance in the naive cascade algorithm that employs no 
subdivision artificially reduces elliptic flow and heats up 
the p± spectra ||, |TT| . 

We modeled central Au+Au collisions at RHIC with 
the minijet initial conditions used in Ref. j^]. The 
evolution started from a longitudinally boost invariant 
Bjorken tube at proper time r 0 = 0.1 fm/c, with locally 
isotropic momentum distribution and uniform pseudora¬ 
pidity g = l/21og[(t + z)/(t — z )] distribution between 
\g\ < 5 with dNg{ro)/dg = 1050. The initial trans¬ 
verse density distribution was proportional to the binary 
collision distribution for two Woods-Saxon distributions, 
while the p± distribution was a thermal fit with T = 700 
MeV as in Ref. ||. 

Pion freezeout results. The freezeout distribution 
d 4 N/d 4 x was defined as the distribution of space-time 
coordinates for the last interaction point of the test parti¬ 
cles. Our strong simplifying assumption is that this point 
is not affected by hadronization. We also neglected reso¬ 
nance contributions to the pion yield. The same 1 g —> l7r 
hadronization model was applied as in Ref. Q motivated 
by parton-hadron duality. 

Fig. 0 shows the pion freezeout distribution d 2 N/rdrdt 
(where r = \/x 2 0 + x 2 s ). Three ranges of the average pair 
transverse momentum are considered. The variation of 
the distributions with transport opacities x = 0.60, 3.01, 
and 7.73 is shown. For dN g /dg = 1050 these correspond 
to a e i = 0.6, 3, and 7.5 mb. Unlike the sharp freezeout 
imposed in hydrodynamical models, the transport th eory 
freezeout is a continuous, evaporation-like process [|l0[ . 
For a given (nonzero) opacity, the larger the p± of the 
particle, the earlier it decouples and the closer it is to the 
surface of the nuclei. Low-p^. particles freeze out from 
the center at late times, while high-pj_ ones escape from 
the surface early. Furthermore, a larger opacity increases 
decoupling times, especially for low-pj_ particles. 

Fig. 0 also demonstrates the importance of using a co¬ 
variant algorithm, such as MPC, for solving Eq. <D- 
For x = 7.73, the covariant freezeout dN/rdrdt distribu¬ 
tion (third row) differs significantly from that obtained 
with the naive noncovariant cascade method, i.e., with¬ 
out particle subdivision (fourth row). At the center, the 
covariant freezeout density decreases with increasing p ±, 
opposite to the monotonic increase shown by the non¬ 
covariant result. Moreover, for the naive algorithm, the 
freezeout density peaks at the center (r « 0) for all p± 
bins shown, while for the covariant result the maximum 
moves out towards the surface as p± increases. 
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FIG. 1: Contour lines for freezeout distributions dN/rdrdt 
from MPC for Au+Au at RHIC as a function of p± and 
transport opacity. Distributions are normalized to unity in 
each p± bin. 


Naturally, Lorentz-violating artifacts also affect the ge¬ 
ometric source radii as shown in Fig. for \ = 7.73. 
Compared to the naive cascade algorithm, the covariant 
one gives smaller Ro and Rl for all I\±. The discrepancy 
is larger for Rl, and for both radii increases as K± de¬ 
creases, leading to about 1 and 2 fm corrections, respec¬ 
tively, at K± fts 0.2 GeV. It is remarkable that for R$ the 
correction changes sign as K± increases, and that con¬ 
sequently the naive cascade method yields Ro/Rs > 1 
for all K l, while the covariant result shows Ro < Rs for 
K± > 2 GeV. In the transport Ro < Rs occurs because 
of strong positive dynamical xo — t correlations. 

Fig. |] shows a comparison to the HBT radii measured 
at RHIC. In the transverse opacity range x ~ 0 - 8 we 
studied, the classical transport results are smaller than 
the observed Ro and Rl- This is in sharp contrast to 
ideal hydrodynamics, which overpredicts both radii 0 . 
The monotonic dependence of Ro and Rl on transport 
opacity suggests that better agreement with data may be 
possible with larger opacity x ~ 20 — 30, which unfortu¬ 
nately are numerically impractical as yet. The need for 
such high opacities is also suggested by the elliptic flow 
V 2 {p±) saturation analysis ||. 

In Fig. [| Rl is obviously small for zero opacity because 
freezeout occurs then at the formation time r = tq. Note 


that R^ ~ r 2 [A(ry — y)] 2 depends on the decoupling time 
and the strength of the g — y correlation. For our ther¬ 
mally correlated initial condition [A(rj — y)} 2 ~ T/m± at 
r 0 = 0.1 fm/c. 

However, as the transport opacity increases, Rl grows 
rapidly because the decoupling time increases as is evi¬ 
dent from Fig. |I]. The largest increase t/t 0 ~ R/t 0 ~ 50 
is for low-p_L partons, which freeze out latest. The ob¬ 
served Rl(K±) is a sensitive probe of the product of 
the freezeout proper time and A (77 — y). Thus a perfect 
inside-outside correlation, i.e., ?/ = y, as assumed in clas¬ 
sical Yang-Mills approaches, cannot be reconciled with 
the RHIC Rl data, without final state interactions. 

In our approach, the main remaining puzzle in Fig. ^|is 
the predicted Rs(K±) ss const « 3 fm that is peculiarly 
independent of the transport opacity and underestimates 
significantly the observed side radius. This suggests that 
Rs is insensitive to the early partonic collective dynam¬ 
ics. The same underestimate of Rg has been found in 
hydrodynamical calculations as welljTj. The Rs problem 
may be related to the assumed longitudinally boost in¬ 
variant dynamics in both approaches. However, it also 
could be related to our neglect of hadronic resonances 
treated in 00 - We also note that for more spherically 
symmetric initial condition, even ideal hydrodynamical 
solutions J 2 l|] exhibit Ro/Rs ~ 1 with larger side radii. 

Conclusions. Using the MPC technique, we investi¬ 
gated the effect of early phase dissipative partonic dy¬ 
namics on the decoupling geometry in heavy-ion colli¬ 
sions in the RHIC energy domain. The pion freezeout 
distribution at midrapidity was found to be sensitive to 
the transport opacity of partons as in |0| . The transport 
freezeout process is similar to evaporation: high-pj_ par¬ 
ticles freeze out early from the surface, while low-p_L ones 
decouple late from the center. For K± < 2 GeV, Ro was 
found to become smaller than Rg indicating that posi¬ 
tive xo — t dynamical correlations are strongly sensitive 
to finite mean free path effects. 

We also demonstrated that the naive cascade algo¬ 
rithms without high particle subdivision lead to large 
numerical artifacts in the freezeout distribution due to 
violation of Lorentz covariance. These artifacts enhance 
the out and long radii, especially at low K±, while in¬ 
crease (reduce) Rs for K _l below(above) « 2 GeV. The 
MPC technique makes it possible to avoid those artifacts. 

While we showed that a decreasing R 0 ut/R S ide < 1 
with K± can arise from covariant parton transport dy¬ 
namics, the momentum scale where this happens is not 
realistic due to the simplified local g —> n hadroniza- 
tion scheme employed as well as the neglect of hadronic 
transport. A consistent explanation of all the differential 
features of HBT correlations will have to include in the 
future a more realistic covariant model of hadronization 
as well as maintain covariance during the hadronic final 
state interactions. 
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FIG. 2: Strong parton subdivision dependence of the pion 
HBT radii as a function of K± for transport opacity x = 7.73. 
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FIG. 3: HBT radii as a function of K± and transport opacity. 
Data from Refs, jl], are shown. 
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